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• Abstract 

' 

q . In view of the successful asymptotic Pade-approximant predictions for higher-loop terms within QCD and 

' massive scalar field theory, we address whether Pade-summations of the MS QCD /3-function for a given number 

of flavours exhibit an infrared-stable fixed point, or alternatively, an infrared attractor of a double valued couplant 
as noted by Kogan and Shifman for the case of supersymmetric gluodynamics. Below an approximant-dependent 
flavour threshold (6 < n/ < 8), we find that Pade-summation j3- functions incorporating [2|1], [1|2], [2 [2], [1 1 3] , 
and [3 1] approximants whose Maclaurin expansions match known higher-than-one-loop contributions to the 
/3-function series always exhibit a positive pole prior to the occurrence of their first positive zero, precluding any 
identification of this first positive zero as an infrared-stable fixed point of the /3- function. This result is shown 
, to be true regardless of the magnitude of the presently-unknown five-loop /3-function contribution explicitly 

appearing within Pade-summation /3-functions incorporating [2|2], [1|3], and [3|1] approximants. Moreover, the 
, pole in question suggests the occurrence of dynamics in which both a strong and an asymptotically-free phase 

^\ . share a common infrared attractor. We briefly discuss the possible relevance of infrared-attractor dynamics to the 

success of recent calculations of the glueball mass spectra in QCD with N c — » oo via supergravity. As n/ increases 
above an approximant-dependent flavour threshold, Pade-summation /3-functions incorporating [2|2], [1|3], and 
[3 1 1] approximants exhibit dynamics controlled by an infrared- stable fixed point over a widening domain of 
the five-loop MS /3-function parameter (/34//3o). Subsequent to the above-mentioned flavour threshold, all 
approximants considered exhibit infrared-stable fixed points that decrease in magnitude with increasing flavour 

1 <""| number. 

a: 

1 Introduction 

Asymptotic Pade-approximant methods have been utilized to estimate higher order contributions to renormalization- 
group (RG) functions within both QCD @, |, || and massive scalar field theory |[ f|, for which such estimates 
compare quite favourably with explicit calculation Q . More recently, such methods have been shown to predict RG- 
accessible coefficients of logarithms within five-loop-order contributions to QCD correlation functions with striking 
accuracy || . These results are all derived from an improvement of Pade-estimated coefficients which incorporates 
the estimated error of Pade-approximants in predicting the n\K~ n n l asymptotic behaviour expected for n th order 
coefficients of a field-theoretic perturbative series [0, ^]. For [AT|M] approximants, such error is seen to decrease 
with increasing N and M jl], g, [tJ] , as well as to favour diagonal and near-diagonal approximants Q . 

Of course, the use of such higher approximants becomes tenable only if the corresponding perturbative series 
is known to sufficiently high order. A known series of the form X^=i specifies all coefficients within [JV|M]- 
approximants to the series only for {TV, M} such that k — N + M; even four-loop calculations (corresponding 
to k = 3 if the leading x 2 behaviour is factored out from the series) serve only to specify [2|1], [1 j 2] , and [0 1 3] 
approximants. Nevertheless, such approximants have been used in conjunction with the anticipated asymptotic 
error to predict the next (five- loop) coefficient i?4 as well as the corresponding diagonal [2\2] approximant to the 
full field-theoretic series. @, |, |, |, § 

The success of such predictions for those cases in which i?4 is known jjl ||, |(| suggests that 
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1. Pade approximants determined from {Ri, Rk} more accurately represent the field-theoretic asymptotic 
series X)j=o Rj x ^ than mere truncation of this series to X^=o Rj&'i 

2. Sufficiently- high Pade-approximants grow arbitrarily close to the function of x represented by the full field- 
theoretic series, as suggested by asymptotic error formulae jjj following from renormalon-estimates of large- j 
coefficients in the series @. 

In the absence of alternatives other than explicit series truncation, such "Pade-summation" || of the full per- 
turbative series may provide a much wanted means for extrapolating such series to the infrared region. We are 
particularly interested in two possible scenarios for infrared dynamics within QCD, either the infrared attractor 
suggested by Kogan and Shifman within the context of supersymmetric gluodynamics ]l0| , or alternatively, dy- 
namics governed by an infrared-stable fixed point. In reference Q], for example, a Pade-summation of the n/ = 3 
QCD /3-function is argued to contain a zero corresponding to an infrared fixed point comparable to that predicted 
by Mattingly and Stevenson 0, [l2|. 

Indeed it is this claim that provides some of the motivation for our present work. In the approach of ref . , 
an infrared-stable fixed point is argued to occur even when rif = 0, in contradiction to a broadening consensus that 
values of n/ even larger than 3 are required for infrared-stable fixed points to occur Jll| |l4|, [l5| |l^]; e.g. nj = 8 
is suggested by a two-loop truncation of the QCD /3-function, with consequences for the phase-structure of QCD 
first explored by Banks and Zaks pq| . 

In the present paper, we utilize the perturbative series for the MS QCD /3-function, now known in full to four- 
loop order |jlQf| , in order to construct [iV|M] Pade summations, which are assumed to provide information about 
the /3-function's first positive zero or pole (this point is further discussed in Section 2). We are able to extend our 
analysis past N + M = 3 by expressing N + M = 4 approximants in terms of the presently-unknown five-loop 
/3 function coefficient, which is treated here as a variable parameter. Among the specific issues we address in the 
sections that follow are: 

1. the existence of a flavour-threshold for dynamics governed by an infrared stable-fixed point, 

2. whether differing Pade-approximants are consistent in predicting infrared properties of QCD, 

3. the dependence of Pade-predictions for /3-function infrared properties on the presently-unknown five-loop 
term, 

4. the size of the infrared fixed point, particularly in comparison with the a* = 7r/4 benchmark value for 
chiral-symmetry breaking |l3|, Q , 

5. the existence of a strong phase of QCD for rif = |l(| and, possibly, for nonzero nj as well, and 

6. the elevation of the true infrared cutoff (mass gap) of nt = 3 QCD to hadronic mass scales (500 - 700 MeV). 

In Section 2, we discuss how Pade-approximants constructed from the known terms of the /3-function series can 
exhibit information about the infrared behaviour of the corresponding couplant. This approach relies upon the 
Pade-approximant remaining closer to the true /3-function than the truncated perturbation series from which the 
approximant is constructed, as discussed above. We conclude Section 2 by obtaining Pade-summation expressions 
for QCD MS /3-functions which incorporate [2|1], [1|2], [2|2], [ 1 1 3] and [3|1] approximants to post-one-loop terms in 
the /3-function series. The latter three approximants are expressed in terms of a variable Ri(= (3^/ ' (3q) characterizing 
the presently unknown 5-loop contribution to the MS /3-function. 

In Section 3, we apply such Pade-summation methods to the /3-function characterizing N c — 3 QCD with 
no fundamental fcrmions. Pade-summation predictions for the infrared structure of rif = QCD in the 't Hooft 
(N c — > oo) limit f20|| are presented separately in an Appendix. For both cases, we find that no Pade-summation 
/3-function supports the existence of an infrared-stable fixed point for the n./ = QCD couplant. Moreover, 
we demonstrate that the infrared behaviour extracted from Pade- summations of the nf = QCD /3-function 
appears to be governed by an apparent /3-function pole, an infrared-attractor of two ultraviolet phases of the 
couplant. This behaviour is in qualitative agreement with that extracted from supersymmetric QCD in the absence 
of fundamental- representation matter fields |U^. We conclude Section 3 with a brief discussion of the possible 
applicability of infrared-attractor dynamics to the glueball spectrum for the N c — > oo case. 
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In Section 4, we extend the analysis of Section 3 to nonzero rif. Specifically, we examine Pade-summation 0- 
functions which incorporate [2|1], [1|2], [2|2], [1|3] and [3 1 1] approximants to post-one-loop terms in the perturbative 
/3-function series. We find that all such approximants exhibit a flavour threshold for the occurrence of infrared 
dynamics characterized by an infrared-stable fixed point. Beneath this threshold, which occurs between 6 and 
9 flavours (depending on the approximant), no infrared-stable fixed point is possible regardless of the magnitude 
of the unknown five-loop term ((34) entering such approximants. Above the threshold, we observe a progressively 
broadening domain of 0a/0o f° r which an infrared-stable fixed point occurs, as well as a decrease in the magnitude 
of such fixed points with increasing rif. 

In Section 5, we focus on the infrared behaviour of the rif — 3 case. We show that [2|2], [1|3] and [3[1] Pade- 
summations of the rif = 3 perturbative /3-function series yield similar infrared dynamics to the "gluodynamic" 
rif = case of Section 3. Such summations are all shown to yield an enhanced mass gap, an infrared boundary to 
the domain of a s somewhat in excess of 500 McV, regardless of the 5-loop contribution to the /3-function series. 
This infrared boundary is shown to be remarkably stable against such 5-loop corrections to the /3-function. 



2 Methodology 

2.1 A Toy f3-Function: 

Pade-approximants to a function whose Maclaurin series is 1 + X)fc=i Rk% k are well known to be valid for a broader 
range of the expansion parameter x than truncations of the series. Consider, for example, the following toy 0- 
function 

~ P A ( X ) = —x 2 [sec{x) — tan(x)] (2.1) 
We have chosen /3a to be asymptotically free, i.e., to have an ultraviolet fixed point at x = 0. Since 



lim (sec(x) — tan(x)) = 



(2.2) 



we have also chosen 0a to have an infrared fixed point at x = ir/2. 0a has a subsequent pole at x = 3ir/2, and 
0a alternates zeros and poles as x increases by subsequent increments of n. The point here, however, is that the 
solution to (2.1) will exhibit the same dynamics as depicted in Fig. 1 for x between zero and n/2, corresponding 
to a freezing-out of the coupling at the x = tt/2 infrared-stable fixed point. 

Suppose, however, that the sum-total of our knowledge of 0a is the first five-terms of this series expansion, 
corresponding to a hypothetical five-loop /3-function calculation: 



0f{x) = -x* 



x 



x 6 5x 4 
2 3~ + ~~2A 



(2.3) 



This truncated series is, of course, not equal to zero at the x = tt/2 infrared fixed point. Rather, when x = w/2, 
each term in the series is seen to be comparable to prior lower-order terms: 



( l\*m 



[1 - 1.571 + 1.233 - 1.292 + 1.269] 



(2.4) 



One would necessarily conclude that the field theoretical calculation leading to (2.3) cannot be extended to large 
enough x to extract information about the infrared properties of 0a{x). 

The series (2.3), however, provides sufficient information to construct a [2|2] approximant to the degree-four 
polynomial within (2.3): 



0T ] (x) = -x* 



-1 X X 

1 2 12 

-1 I X X 2 

1 "r" 2 12 



(2.5) 



Equation (2.5) is obtained by requiring that the degree-2 numerator and denominator polynomials of the [2|2] 
approximant be chosen so as to yield a Maclaurin expansion whose first five terms reproduce the five terms in (2.3). 
One can easily verify that /3^' 2 ' remains closer to 0a (2.1) over a much larger range of x than A 5 \ as given in 
(2.3). This range is inclusive of the first zero of 0a- 0a^ h as a positive zero at x = 1.583, quite close to 0a's true 
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zero at 7r/2. Moreover, the denominator in (2.5) remains positive over the entire range < x < 1.583, guaranteeing 
that x — 1.583 is infrared-stable (a sign change would render this fixed-point ultraviolet-stable). Thus, Pade- 
improvement of the information in (2.3) provides a means for extracting information about the infrared properties 
of Pa that is otherwise inaccessible from the "five-loop" expression. 



2.2 (3-Function Poles 

It should also be noted that (2.5) predicts that a pole at x = 7.58 follows the zero at 1.583 without a second 
intervening zero. This result is qualitatively similar to the true behaviour of (3a, which acquires a pole at x = 3ir/2 
subsequent to the zero at n/2 without any additional intervening zeros. However, accuracy in predicting this pole, 
as well as any subsequent zeros or poles, is clearly beyond the scope of (2.5), the [2|2] Pade-summation of (3a- 

We have seen, however, that Pade methods do provide a window for viewing leading /?-function singularities that 
would otherwise be inaccessible. One cannot automatically dismiss the possibility of such singularities occurring 
within QCD /3-functions. For example, the /3-function of SU(N C ) SUSY gluodynamics, which is known exactly if 
no matter fields are present, exhibits precisely such a zero Ell: 



(3{x) 



W c x 2 



1 - N c x/2 



7T 



(2.6) 



Eq. (2.6), which can be derived via imposition of the Adler-Bardeen theorem upon the supermultiplet of the 
anomalies [22 , implies the existence of a strong ultraviolet phase (the upper branch of Fig. 2) when the couplant 
x is greater than the /3-function pole at 2/N c JhJ. Interestingly, the /3-function (2.6) is itself a [0|1] approximant 
once the leading — 3N c x 2 /A coefficient is factored out. 

To demonstrate how Pade summation provides a window for extracting possible pole singularities in true (3- 
functions, we consider a second toy example 



dx 



(3b{x) = x 2 [sec(x) + tan(x)] 



(2.7) 



(3b(x) is asymptotically free, but has a positive pole at x = n/2 prior to its first zero at x = 3tt/2. This zero 
is an ultraviolet stable fixed point because of the overall sign change associated with passing through the pole at 
x = 7r/2. The behaviour of x(/i) for x < 3n/2 is schematically depicted in Fig 2, with x — 7r/2 corresponding to 
He the minimum allowed value of \i (assuming the couplant x is real). 

Such infrared structure is not at all evident in the "five-loop" approximation to (2.7) 



(3%\x) 



1 



x x 5x 
~2 + T + ~2A 



(2.8) 



an expression which ceases to be close to the true (3- function (2.7) for values of x substantially smaller than x = n/2. 
However, one can obtain a [2|2] approximant directly from the truncated series in (2.8) 



/4 2|2I (*) 



12-1 



12-1 



(2.9) 



whose Maclaurin expansion yields (2.8) for its first five terms. The first denominator zero of (2.9) is at x = 1.583, 
in good agreement with the positive pole of (2.7) at x — tt/2. Moreover, the first denominator zero of (2.9) precedes 
all (positive) numerator zeros, thereby eliminating the possibility of the pole being preceded by an infrared fixed 
point. Thus, (2.9) and the true /3-function (2.7) predict very similar dynamics between the ultraviolet-stable fixed 
point at x — and the infrared-attractor pole at x — tt/2. By contrast, the "five-loop" /3-function (2.8) can only 
reproduce the true running couplant in the ultraviolet region where x is near zero. 



2.3 Pade- Improvement and Infrared Behaviour 

It is to be emphasized that the examples presented above demonstrate how Pade-improvement may provide infor- 
mation about the infrared region that a truncated perturbative series cannot. There is no way, of course, to prove 
that the first positive zero or pole of a given Pade-summation is indeed the first zero or pole of the true /3-function. 
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In the absence of methodological alternatives, however, we will explore below the consequences of assuming this to 
be the case. Corroboration of such an assumption relies ultimately on an explicit comparison of next-order terms 
calculated for the MS QCD /3-function series, to Pade-predictions for these terms {e.g., ref ||). There is reason to 
be encouraged, however, by the success already demonstrated for Pade predictions of the known five-loop content 
of the massive scalar field theory /3-function |[| [|. Similar successes are obtained in predicting RG-accessible 
coefficients within the five loop contributions to QCD vector and scalar fermionic-current correlation functions, as 
well as within four- loop contributions to the scalar gluonic-current correlation function |5| . 

The general approach we take is to replace a fc-loop truncation of the asymptotic /3-function series 

(3^ k \x) = -p x 2 (1 + Rix + ...Rk^- 1 ) (2.10) 

with an expression incorporating the corresponding [iV|JVf] Pade- approximant {N + M = k — 1): 



l*\M] M _ ( l + aix + ... + a N x 



N 



The N+M coefficients {ai, ajy, bi, &m} are completely determined by the requirement that the first k terms in 
the Maclaurin expansion of (2.11) replicate (3^{x) (2.10). We then examine /3l JV 'l Af ] in order to determine whether 
or not it is supportive of an infrared-stable fixed point, as in Fig. 1, or an infrared-attractor pole, as in Fig. 2. 

If the first positive zero of the degree-N polynomial in the numerator of (2.11) precedes any positive zeros of the 
degree-M polynomial in the denominator, that first numerator zero corresponds to the infrared-stable fixed point 
at which the couplant x freezes out in Fig. 1. Alternatively, if the first positive zero in the denominator of (2.11) 
precedes any positive zeros in the degree-M numerator polynomial, that first denominator zero corresponds to the 
infrared-attractor pole common to both couplant phases of Fig. 2. 

2.4 A Pade Roadmap 

It will prove useful to tabulate those formulae required to obtain Pade approximants (2.11) whose Maclaurin 
expansions reproduce the truncated series (2.10) for the MS /3-function. Values for (3q, R%, R2 and R3 for the 
/3-function, as defined by (2.10), are tabulated in Table |l|. Corresponding [2|1] and [1|2] approximants to the 
truncated series 1 + R\x + R2X 2 + R3X 3 within (2.10) are given by the following formulae: 



n f 


A) 


Ri 


R2 


R 3 





11/4 


51/22 


2857/352 


41.5383 


1 


31/12 


67/31 


62365/8928 


34.3295 


2 


29/12 


115/58 


48241/8352 


27.4505 


3 


9/4 


16/9 


3863/864 


20.9902 


4 


25/12 


77/50 


21943/7200 


15.0660 


5 


23/12 


29/23 


9769/6624 


9.83592 


6 


7/4 


13/14 


-65/224 


5.51849 


7 


19/12 


10/19 


-12629/5472 


2.42409 


8 


17/12 


1/34 


-22853/4896 


1.00918 


9 


5/4 


-3/5 


-1201/160 


1.97366 


10 


13/12 


-37/26 


-41351/3744 


6.44815 


11 


11/12 


-28/11 


-49625/3168 


16.3855 


12 


3/4 


-25/6 


-6361/288 


35.4746 


13 


7/12 


-47/7 


-64223/2016 


71.6199 


14 


5/12 


-113/10 


-70547/1440 


145.373 


15 


1/4 


-22 


-2823/32 


332.091 


1G 


1/12 


-151/2 


-81245/288 


1309.98 



Table 1: The known coefficients /3 and R\ 2 3 appearing in the QCD /3-function (2.10) for nj — — 16, as calculated 
in M. 
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f}™(x) = -0 o x 



1 + a\x + a,2X 2 



1 + bix 
a>i — (R1R2 — Rz)/R2, 
a-2 = (R 2 — RiRa)/Ri, 
b\ = —R3/ R2, 



and 



f3m( x ) = -p x 2 



1 + a\x 



1 + + b2X 2 

h = (R 3 - R l R 2 )/{Rl - R 2 ), 
b 2 = (Rl - RiR 3 )/(Rl - R2), 
ai = Ri + bi. 



(2.12a) 

(2.12b) 
(2.12c) 
(2.12d) 

(2.13a) 

(2.13b) 
(2.13c) 
(2.13d) 



We do not consider the [0|3]-approximant, as this approximant has no possible numerator zeros and therefore 
cannot lead to an infrared fixed point. The [3 10] -approximant is, of course, the truncated series itself. The 
"diagonal-straddling" [1|2] and [2|1] approximants are the only N + M = 3 approximants for which both infrared- 
stable fixed points (Fig. 1) and infrared-attractor poles (Fig. 2) are possible, depending on the specific ordering of 
positive numerator and denominator zeros in (2.12a) and (2.13a). Hence, it is these approximants we will study in 
subsequent sections. 

i?4, the "next-order" coefficient of the QCD /3-function (2.10), is not presently known. Nevertheless, one 
can construct "diagonal-straddling" [2|2], [1|3] and [3|1] approximants with R4 taken to be an arbitrary pa- 
rameter, by utilizing the values for {R±, R2, R3} given in Table |l|. One finds for the truncated series j3^\x) 
= —[3qx 2 [l + R\x + R2X 2 + R3X 3 + R4X 4 ] the following Pade-approximant /3-functions: 



0™{x) 
h = 



-Pox 2 



1 + CL\X + 020; 



1 + b\X + 

(R1R4 — R2R3) I \R\ — R1R3) 



^2 = {R3 — R2Ra)/{R2 — R1R3), 

a\ = Ri + 61, 
a 2 = R\b\ + b 2 + i?2; 



1 + d\X 



1 + b\x + 62X 2 + 63a; 3 
(i? 2 i?2 — R\ — R1R3 + R4) 



&3 



{2R 1 R 2 -Rf- R 3 ) ' 

(i? 2 i?3 — i?ii?2 + R2R3 ~ R4R1) 

~ 2R X R 2 -RI-R3 ' 

[R 2 + i?f — 2R1R2R3 + Ri(R\ — R2)} 
2R X R 2 -Rf-R 3 

ai = Ri + bi. 



(3 m (x) = -p x 



1 + d\X + CL2X + (I3X 3 



1 + hx 
a% = -Ri — R4/ i?3, 
o-2 = R2 — RiRi/ R31 

a 3 = R3 — R2R4/R31 



2.14a) 

2.14b) 
2.14c) 
2.14d) 
2.14e) 

2.15a) 

2.15b) 

2.15c) 

2.15d) 
2.15e) 

2.16a) 

2.16b) 
2.16c) 
2.16d) 
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h = -R4/R3. 



(2.16e) 



Given the known values of {Ri, R2, R3} tabulated in Table [I], we have tabulated /jl 2 ' 2 ! 's coefficients {01, 02, &i, 62} 
in Table |^ for all n/-values for which x — is an ultraviolet-stable fixed point. These coefficients are all linear in the 
unknown parameter R4. A similar tabulation of the coefficients within /3W 3 ] and Z^ 3 ' 1 ! is presented in Tables || and 
[|, respectively. Because all such Pade-coefficients are i?4-dependent for a given choice of n/ , one might expect to 
find infrared stable fixed point behaviour (Fig. 1) for some range of R4, infrared-attractor pole behaviour (Fig. 2) 
for another range of Ri, and (possibly) some regimes of Ri for which there are neither numerator nor denominator 
zeros. Surprisingly, we find in Section 4 that infrared-stable fixed point behaviour as in Fig. 1 does not occur for 
any of the diagonal-straddling approximants (2.12-15) until nf reaches a threshold value, regardless of Ri. This 
threshold depends on the approximant considered, but is greater than or equal to 6 for all [iV|Af] approximants 
discussed above. 

n f ax a 2 61 b 2 



13.4026-0.0762155 Ri 

1 116019-0.085085&R4 

2 9.50934-0.0941220i? 4 

3 7.19456-0.102610i? 4 

4 4.84008-0. 110684i? 4 

5 2.67929-0.123291^4 

6 0.610851-0.1842367? 4 

7 1.907466+0.129932i? 4 

8 0.245912+0.001351797? 4 

9 -0.342477-0.0104297^4 

10 -0.880095-0.0108500i?4 

11 -1.65139 - 0.00886659J? 4 

12 -2.93401-0.006555097?4 

13 -5.18889-0.004488997?4 

14 -9.53837-0.002795077?4 

15 -20.0584-0.001458067? 4 

16 -73.4295-0.000423006i? 4 



-22.9153+0.0901663i? 4 
-19.0068+0.0911037i? 4 
-15.0709+0.0875660i? 4 
-11.3292+0.07564387?4 
-8.18415+0.04858877? 4 
-6.19674-0.0112453i? 4 
-6.62748-0.2286507?4 
-0.130346+0.638145i? 4 
-4.61451+0. 21457LR 4 
-7.59305+0.136738i? 4 
-11.5003+0.099648^4 
-17.0050+0.0771335i?4 
-25.2430+0.0620604i? 4 
-38.6692+0.05143897? 4 
-63.6700+0.04370237?4 
-123.626+0.0379240i? 4 
-428.807+0.0335175i? 4 



-0.0762155i? 4 + 11.0844 
-0.0850858i? 4 +9.44058 
-0.0941220i? 4 +7.52659 
-0.102610i? 4 + 5.41678 
-0.110684i? 4 +3.30008 
-0.123291i? 4 +1.41842 
-0.184236i? 4 -0.317721 
0.12993187? 4 +1.38115 
0.00135179i? 4 +0.216500 
-0.0104297i? 4 +0.257228 
-0.0184997.R4 +0.542982 
-0.00886659i? 4 + 0.894061 
-0.00655509 J? 4 +1.23265 
-0.00448899i? 4 +1.52540 
-0.00279507i? 4 +1.76163 
-0.00145806J? 4 +1.94165 
-0.000423006i? 4 + 2.07047 



0.266848i? 4 - 56.7275 
0.274999i? 4 -46.3959 
0.274187i? 4 -35.7703 
0.258062i? 4 -25.4301 
0.219042i?4-16.3139 
0.144208i? 4 - 9.45997 
-0.0575739i? 4 -6.04228 
0.569760i? 4 +1.45066 
0.214531i? 4 +0.0468084 
0.1304807?4+0.0677118 
0.0842074i? 4 +0.317008 
0.0545640i? 4 + 0.935218 
0.0347475i?4+1.97982 
0.0212985i? 4 +3.42939 
0.01211807?4+5.22737 
0.00584673i? 4 +7.30916 
0.00158053i? 4 +9.61460 



Table 2: Coefficients ai, 2 and 61,2 of (3^ 2 \ the [2|2] Pade-approximant (2.13) to the QCD (3 function for nj = — 16. 
The coefficients are all linear in Ri, the (presently-) unknown five- loop term in (2.10). 



3 Gluodynamics 

In this section we consider conventional (N c = 3) QCD with nf = 0. The N c — > 00 case is considered separately 
in an Appendix. This "gluodynamic" limit is of particular interest as a possible projection (without gluinos) of 
supersymmctric QCD in the absence of fundamental-representation matter fields (rif — SQCD), a theory for 
which the /3-function is known to all orders of perturbation theory pjj. The nj — SQCD /3-function (2.6) does 
not exhibit an infrared-stable fixed point; rather, it exhibits the dynamics of Fig. 2 in which a /3-function pole 
serves as an infrared-attractor, both for a weak asymptotically-free phase, as well as for a strong phase of the 
now double- valued couplant JlO[ ] . Such dynamics differ fundamentally from those of Fig. 1 anticipated from an 
infrared-stable fixed point, which has been argued elsewhere [ [l2"[ to occur for QCD even in the nj — gluodynamic 
limit. Pade-approximant estimates of higher order terms in QCD /3-functions j[) and correlators have in fact 
proven to be most accurate in the n/ = case in which "quadratic-Casimir" effects are minimal jl|. Thus, the 
methodological machinery of the previous section may be particularly well-suited to shed insight on whether the 
dynamics of Fig. 1 or Fig. 2 characterize QCD's gluodynamic limit. 
For n f = 0, the four-loop ~MS QCD /3-function is H) 

(3 {A) {x) = -yi 2 [l + 2.31818a; + 8.11648a; 2 + 41.5383a; 3 ] , (3.1) 
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n f 


<2l 


01 


02 


03 





9.56239-0.061 1053/14 


7. 24421-0. 061 1053/i4 


-24.9099+0. 141653/C4 


-42.5901+0. 167582it4 


1 


8. 51104-0. 0702708/14 


6.34975-0.0702708/14 


-20.7090+0. 151876it4 


-33.9265+0. 162617i<4 


2 


7.25658-0. 0810329it4 


5. 27382-0. 0810329it4 


-16.2327+0. 160669it4 


-25.7265+0. 149477Jt4 


3 


5. 80845-0. 0933552/t4 


4. 03067-0. 0933552it4 


-11.6367+0. 165965it4 


-18.3242+0.122349it4 


4 


4.24716-0. 107164it4 


2.70716-0. 107164it4 


-7.21667+0. 165032it4 


-12.2028+0.07244688it4 


5 


2.76704-0. 123131it4 


1.50617-0. 123131/u 


-3.37387+0. 155253it4 


-7. 80319-0. 0141605it4 


6 


1.72453-0. 145814/i4 


0. 795958-0. 145814/t4 


-0.448926+0. 135399i<4 


-4.87066-0. 168040it4 


7 


1 .97486-0. 200029it4 


1. 44855-0. 200029-K4 


i r a r r a i r\ 1 at ato t~) 

1 .54554+0. 105278it4 


0. 105614-0. 517062it4 


8 


17.0270-0. 778954i?4 


16. 9976-0. 7789537917J<4 


4. 16776+0. 0229104it4 


78. 2076-3. 63659it4 


9 


-8.581 12+0. 137934i?4 


-7.981 12+0. 137934i?4 


2. 71758+0. 08276047?4 


-60. 2514+1. 08502i?4 


10 


-6. 27351+0. 03588297?4 


-4. 85044+0. 0358829i?4 


4. 14206+0. 0510641i?4 


-54. 12482+0. 468980i?4 


11 


-6.36697+0.01252297? 4 


-3.82151+0.0125229i? 4 


5.93697+0.03187657? 4 


-61.1352+0.277305i?4 


12 


-7.44146+0.004526527? 4 


-3.27479+0.004526527? 4 


8.44183+0.0188605.R4 


-72.6300+0.178562J?4 


13 


-9.70448+0.00151777^4 


-2.99019+0.00151777i? 4 


11.7796+0.010191^4 


-87.7855+0. 116775.R4 


14 


-14.2164+0.000415849i? 4 


-2.91637+0.000415849i? 4 


16.0360+0.004699097? 4 


-107.043+0.0734726% 


15 


-25.0410+7.04348- lO -5 ^ 


-3.04098+7.04348- 10" 5 i? 4 


21.3172+0.00154956i? 4 


-131.384+0.04030417? 4 


16 


-78.8684+2.12019-10- 6 i? 4 


-3.36839+2.12019-10" 6 i? 4 


27.7873+0.000160074i? 4 


-162.269+0.01268377? 4 



Table 3: Coefficients a\ and 61,2,3 of /?^' 3 l, the [1 1 3] Pade-approximant (2.14) to the QCD /3-function for n/ = — 16. 
i?4 is the unknown five-loop contribution to the /3-function. 



as is evident from substitution of the nf — Table gentries into (2.10). This /^-function is sufficient to determine 
the [1 1 2]— and [2|l]-approximant /^-functions via (2.12) and (2.13): 



1 - 2.7996x - 3.7475a; 
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P 



11 2 



1 - 5.1178a; 
1 - 5.9672a; 



1 - 8.2854a; + 11.091a; 2 



(3.2) 



(3-3) 



In both (3.2) and (3.3), the first positive denominator zero (%d) precedes the first positive numerator zero (x n ). 
We find from (3.2) that Xd = 0.195 < x n = 0.264, and from (3.3) that Xd = 0.151 < x n = 0.168. As discussed in 
Section 2, this < Xd < x n ordering of positive zeros is consistent with dynamics in which Xd serves as in infrared 
attractor for both a strong and a weak asymptotically- free phase (Fig. 2). The first positive numerator zero x n , if 
taken seriously, necessarily corresponds to an ultraviolet-stable fixed point because of the /3-function sign-change 
occurring as x passes through Xd- 

To test the stability of these conclusions against higher-than-four loop corrections, we add an arbitrary "five- 
loop" correction to (3.1): 

P {5 \x) = /3 (4) (x) - ^-x 2 [R A x 4 }. (3.4) 

The five-loop /3-function (3.4) determines [2|2L [1|3], and [3|1] Pade-approximant /3-functions via (2.14), (2.15) and 
(2.16). These can be read off Tables || and we list them explicitly here to facilitate the analysis which follows: 



f3^ 2 \x) = ^x 2 



1 + (13.403 - 0.076216^4)2; - (22.915 - 0.090166^ 4 )a; 2 
1 + (11.084 - 0.076216i? 4 )a; + (56.728 - 0.26685i? 4 )a; 2 



(3.5) 



/3 [l|3] (a; ) = _^ a .2 



1 + (9.5624 - 0.061 105i? 4 )x 



1 + (7.2442 - 0.061105i? 4 )x - (24.910 - 0.14165 J R 4 )x 2 - (42.590 - 0.16758i? 4 )a; 3 



(3.6) 



0^(x) 
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1 + (2.31818 - 0.024074i? 4 )x + (8.1165 - 0.055808i? 4 )x 2 + (41.538 - 0.19540i? 4 )a; 



3 1 



1 - 0.024074^ 4 a; 



(3.7) 



Figure 3 exhibits a plot of x n , the first positive zero of (3.5), and Xd (the first positive denominator zero of 
(3.5), as a function of the independent variable i? 4 . Such positive zeros are seen to occur for all R4. Moreover, 
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«1 



a 2 



«3 



2.31818-0.0240742i? 4 

1 2.16129-0.0291294i? 4 

2 1.98276-0.0364292i? 4 

3 1.77778-0.0476412724 

4 1.54000-0.0663747i? 4 

5 1.26087-0.101668i? 4 

6 0.928571-0. 181209 R 4 

7 0.526316-0.412526i? 4 

8 0.0294118-0.990906i? 4 

9 -0.600000-0.506673J? 4 

10 -1. 42308-0. 155083i? 4 

11 -2.54545-0.0610294i? 4 

12 -4.16667-0.0281892i? 4 

13 -6.71429-0.0139626i? 4 

14 -11.3000-0.006878847?4 

15 -22.0000-0.003011227?4 

16 -75.500-0.000763368i? 4 



8.11648-0.0558083i? 4 
6.98533-0.0629572i? 4 
5.77598-0.07223027? 4 
4.47106-0.08469547? 4 
3.04764-0. 102217i? 4 
1.47479-0. 128190i? 4 
-0.290179-0.168265i? 4 
-2.30793-0.2171197?4 
-4.66769-0.0291443i? 4 
-7.50625+0.304004i? 4 
-11.0446+0.220695i? 4 
-15.6645+0. 155348i? 4 
-22.0868+0. 117455i? 4 
-31.8566+0.09374897?4 
-48.9910+0.077730i? 4 
88.2188 + 0.0662468^4 
-282.101+0.0576343i? 4 



41.5383-0.195397J? 4 
34.3295-0.2034797? 4 
27.4505-0.2104147? 4 
20.9902-0.2130077?4 
15.0650-0.202286i?4 
9.83592-0. 149939.R4 
5.51849+0.0525830J? 4 
2.42409+0.952081J? 4 
1.00918+4.62524i? 4 
1.97366+3.803227?4 
6.44815 + 1.71283i? 4 
16.3855+0.9559937? 4 
35.4746+0.622609i? 4 
71.6199+0.4448017? 4 
145.373+0.337001i? 4 
332.091+0.2656467?4 
1309.98 + 0.215347i? 4 



-0.0240742i? 4 
-0.02912947? 4 
-0.0364292i? 4 
-0.0476412i? 4 
-0.06637477? 4 
-0.101668i? 4 
-0.181209i? 4 
-0.412526i? 4 
-0.990906i? 4 
-0.506673^4 
-0.1550837? 4 
-0.0610294i? 4 
-0.0281892i? 4 
-0.0139626i? 4 
-0.006878847?4 
-0.00301122i? 4 
-0.0007633687? 4 



Table 4: Coefficients 01,2,3 and b\ of Z^ 3 ' 1 ', the [3 1 1 ] Pade-approximant (2.15) to the QCD /3-function for n/ = 0- 
i?4 is the unknown five-loop contribution to the /3-function. 



16. 



the first positive denominator zero is seen to precede the first positive numerator zero over the entire range of R4. 
This last result confirms that the Fig. 2 dynamics predicted from [1|2] and [2|1] approximants appear to be stable 
against 5-loop corrections of arbitrary magnitude. 

These results are corroborated by and /3[ 3 W. Figure 4 plots the first positive numerator and denominator 
zeros of (3.6) against R4 for (3^ 3 \ as given by (3.6). A positive numerator zero (x n ) exists only if i?4 > 150, 
whereas at least one positive denominator zero (xd) occurs for all R4. Once again, however, Xd precedes x n over 
the entire range of R4, suggesting that Xd exists as an infrared-attractor for all values of R4. Moreover, the ordering 
< Xd < x n , over all values of i?4 for which x n exists, precludes any identification of x n with an infrared-stable 
fixed point. 

Figure 5 plots x n and Xd against i?4 for /3' 3 ' 1 !, as given in (3.7). This case is perhaps the most interesting 
of all, as a positive denominator root Xd is possible only if R4 > 0, as is evident from the denominator of (3.7). 
Figure 5 shows that Xd continues to precede x n for all positive values of R4. Moreover, when R4 is negative (and 
a positive pole Xd is no longer possible), one sees from (3.7) that the numerator polynomial coefficients are all 
positive-definite, precluding any possibility of positive numerator roots. We thus see that an infrared-stable fixed 
point for /jl 3 ! 1 ! is unattainable, even for the R4 < region for which no denominator zero occurs at all. 

Thus, there does not exist any valid N + M = 4 approximant (M > 0) to the nf = QCD /3-function which 
supports an infrared-stable fixed point, regardless of the magnitude of the presently unknown 5-loop term entering 
such approximants. [] Moreover, the Fig. 2 dynamics following from a positive /3-function pole preceding any 
/3-function zeros appear to be upheld not only by [2|1]- and [l|2]-approximant /3-functions, but also by [2|2]-, [1|3]- 
and (when R4 > 0) [3|l]-approximant /3-functions as well. 

One can also utilize asymptotic Pade approximant methods to estimate the magnitude of the infrared attractor 
(Fig. 2), corresponding to the first positive denominator zero of all the approximants considered so far. An 
asymptotic error formula Q enables one to obtain an asymptotic Pade-approximant prediction (APAP) of R4 from 
the three preceding terms {i?i, R2, R3} in the /3-function series H: n 



R 



APAP 



Rl [Rl 



R1R2R3 — 2i? 3 i? 3 ] 



i?2 [2-Ro — R\R3 ~ ^1^2] 



(3.8) 



Utilizing the values {R 1 ,R 2 ,R 3 } listed in Table |l] for n f = 0, one finds from (3.8) that Rf PAP = 302.2. If one 

1 /3l°l 4 l cannot have a nonzero fixed point, as this approximant has no positive numerator zeros at all. The M = case [/3 I °] ] , 
corresponding to the truncated series itself, exhibits a positive zero identifiable with an infrared fixed point only if R4 < 0, in which 
case the negative five-loop term must be equal in magnitude to the sum of the preceding (positive) one-through-four-loop terms. 

2 This formula summarizes the content of the algorithm presented in Sections 2 and 5 of ref. [nj 
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inserts this value of R4 into the [2 1 2] - , [1 1 3] - and [3|l]-approximant /3-functions (3.5), (3.6) and (3.7), one finds 
reasonable agreement from all three (3- functions as to the location of the first positive denominator zero, i.e. the 
infrared-attractor of Fig. 2. This zero is seen to occur at Xd = 0.11 for fft 2 ^ and /?P-I 3 1, and at x c i = 0.14 for /jl 3 ' 1 !, 
corresponding to values of a s between 0.35 and 0.44. A similar i?4-estimate (R4 = 270) can be obtained for the 
rif = case via the weighted asymptotic Pade approximant procedure (WAPAP) delineated in Section 5 of ref. 
This procedure expresses the coefficient P4 = P0R4 as a degree-4 polynomial in n/, as would be obtained from an 
explicit perturbative calculation. ^ We reiterate that the rif = case minimizes unknown quadratic Casimir effects 
and avoids entirely the large uncertainties associated with the cancellation of large n/-dependent terms within /Jj 
(as evident from eq. (5.5) of Jl|) that characterize rif > estimates based on APAP and WAPAP methods. 

We conclude this section by reiterating that the ordering < Xd < x n of positive denominator and numerator 
zeros of n/ = Pade-approximant /3-functions suggests the existence of a double- valued couplant (Fig. 2), as 
already seen in SUSY gluodynamics [|l0). Such a scenario is seen to decouple the infrared region (^ < fi c ) from 
the domain of purely-perturbative QCD, the domain of (real) a s . Such a scenario is perhaps also indicative of an 
additional phase distinguished by strong coupling dynamics at short distances ^] Taking the Pade-approximants of 
the /3-function seriously, one concludes that there is an ultraviolet-stable fixed point in that phase. Recall that in 
SUSY gluodynamics, such a fixed point is at x = 00 [ jTo| . Of course, it is impossible to assert that this ultraviolet- 
stable fixed point survives in the exact /3-function of (non-SUSY) gluodynamics. For example, it may be replaced 
by an ultraviolet Landau pole. 

Regardless of these considerations, the picture with an infrared attractor seems plausible and self-consistent. In 
such dynamics, both phases may share common infrared properties | jTo| . The presence of two phases has implications 
meriting further exploration. Such dynamics are shown in the Appendix to characterize QCD in the N c — > 00 limit 
for all but the aforementioned [3|l]-case with i? 4 negative. Indeed, such dynamics may prove pertinent to the 
unexpected agreement between the glueball mass spectra obtained via lattice methods 1 23 and those obtained via 
supergravity wave equations in a black hole geometry [Q following from conjectured duality to large-N c gauge 
theories p5| . an agreement obtained despite the large bare coupling constant necessarily utilized in the latter 
approach. 

In the section which follows, we extend the analysis of the QCD /3-function to nonzero rif values. We specifically 
seek to address whether Pade-methods indicate a flavour-threshold for infrared-stable fixed points. However, we 
also seek insight as to whether there is any evidence for a strong phase of QCD when n/ > 0, as such a phase could 
(conceivably) provide a dynamical mechanism for electroweak symmetry breaking. 



4 QCD with Fermions 

In this Section, we repeat the analysis of the previous section with rif < 16 fermion flavours. The rif = 16 case 
is the maximum number of flavours consistent with asymptotic freedom: a /3-function whose sign is negative as 
x 0+ 



4.1 Four-Loop-Level Results 

We first consider the /3-function (2.12) incorporating a [2 1 1 ] approximant to describe post-one-loop behaviour. This 
/3-function is fully determined by the known two-, three-, and four-loop terms {Hi, R2, R3} tabulated in Table [l] for 
??/ = {1, 2, 16}. The values of the coefficients {ai, ai, b\] characterizing /3I 2 ' 1 ] in (2.12a) are tabulated in Table |[ 
Also tabulated in the table are values for x n , the first positive zero of the numerator 1 + a\x + 02a; 2 , as well as Xd, 
the zero of the denominator 1 + b\x. Blank entries for x n , Xd correspond to cases where no positive zero exists. We 
see from the table that a positive denominator zero precedes the first positive numerator zero for rif < 5, in which 
case that denominator zero serves as an infrared attractor for both a strong and an asymptotically-free ultraviolet 
phase of the QCD couplant. The denominator zero becomes negative for rif > 6; nevertheless, an infrared-stable 
fixed point (associated with a positive numerator zero) does not occur until rif — 7, as no positive numerator zero 
exists for rif — 6. The value of x n — a s /Tt associated with this infrared-stable fixed point decreases as rif increases, 

3 Coefficients (3^ listed in |Q must be divided by 4 fc+1 to correspond to our normalization of /3-function coefficients: our f3(x) = 
dx/d(log (i 2 ) with x = a/71-. 

4 lt should be noted here that the presence of an infrared-stable fixed point a* also implies a possible strong phase of the couplant 
, in addition to the asymptotically-free phase exhibited in Fig. 1. If a([i) > ot*[(l > 0], then a approaches a* from above as fi 
approaches zero from above. 
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as anticipated in other work [^3| [Ij]. If we regard a cr = 7r/4 (x n = 1/4) as the threshold value for chiral symmetry 
breaking p6[ , we see that the conformal window for QCD is predicted to begin at rif — 11. For rif = {7 — 10}, the 
infrared-stable fixed point x n is not expected to govern infrared dynamics. Chiral-symmetry breaking should occur 
before the couplant reaches x n , and the (now-massive) fermions are expected to decouple from further infrared 
evolution of a s @ 0. 



/ I f 


ai 


a-2 


6i 




Xd 





-2.79959 


-3.74745 


-5.11777 


0.26394 


0.19540 


1 


-2.75323 


-3.63638 


-4.91452 


0.26820 


0.20350 


2 


-2.76977 


-3.64714 


-4.75253 


0.26710 


0.21041 


3 


-2.91691 


-3.87504 


-4.69468 


0.25586 


0.21301 


4 


-3.40349 


-4.56534 


-4.94349 


0.22557 


0.20229 


5 


-5.40850 


-6.93442 


-6.66937 


0.15435 


0.14994 


6 


19.9461 


17.3690 


19.01757 






7 


1.57664 


-1.75513 


1.05033 


1.3275 




8 


0.245617 


-4.66133 


0.216205 


0.49027 




9 


-0.337065 


-7.66401 


0.262935 


0.33990 




10 


-0.839249 


-11.8754 


0.583828 


0.25699 




11 


-1.49942 


-18.3271 


1.04603 


0.19624 




12 


-2.56052 


-28.7791 


1.60615 


0.14716 




13 


-4.46609 


-46.9517 


2.24819 


0.10593 




14 


-8.33265 


-82.5220 


2.96735 


0.070620 




15 


-18.2356 


-171.036 


3.76441 


0.039903 




16 


-70.8563 


-632.698 


4.64367 


0.012678 





Table 5: The coefficients a±, a-i and b\ in /3' 2 ' 1 ', as given in Eq. (2.12) for rif = 0— 16. The column x n tabulates the 
first positive zero of the numerator. The column Xd tabulates the denominator zero, if positive. For those values 
of rif for which real positive numerator (denominator) zeros do not exist, the corresponding entries for x n (xd) are 
left blank. The table shows < Xd < x n for rif < 5, corresponding to Fig. 2 type dynamics for this range of rif. 
When nf > 7, x n corresponds to an infrared-stable fixed point that decreases as n/ increases. 

Qualitatively similar conclusions arc obtained from the analysis of /jW 2 ! [eq. (2.13)], although the corresponding 
values of n/-thresholds for various infrared properties differ somewhat from those of the [2|l]-approximant case. In 
Table ||the values of the constants {a\, b%, 62} characterizing /jl 1 ! 2 ! in (2.13a) are tabulated using the known values 
for {i?i, i?2, R3} listed in Table [j]. When positive, the zero of the numerator 1 + a\x fails to precede positive zeros 
of the denominator 1 + 61a; + b-ix 2 until rif = 9 (positive denominator zeros cease occurring after rif = 6, but the 
numerator zero is negative when 5 < rif < 8). Consequently, rif = 9 is the flavour-threshold for identification of 
x n as an infrared-stable fixed point. As before, this fixed point decreases with increasing rif. Its magnitude does 
not fall below the x n = 1/4 threshold for chiral-symmetry breakdown until rif — 12, corresponding to a previous 
prediction jl3) of the threshold for QCD's conformal window, consistent with the qualitative picture presented in 
ref. jlJJ . Specific predictions for x n from Tables ^ and ^ agree quantitatively only within this 12 < rif < 16 window. 
Although the intermediate range of rif for which x n > 1/4 differs between the two approximants (7 < rif < 10 
for Z?! 2 ' 1 !; 9 < rif < 11 for Z^ 1 ' 2 !), it is nevertheless significant that such a range exists for both cases. In Table 
H a positive denominator zero is seen to precede any positive numerator zeros when rif < 6. This behaviour 
corresponds to the infrared dynamics suggested by Figure 2. Once again, qualitative agreement is seen to occur 
between the [2|1]- and [l|2]-approximant (3- functions insofar as both 0- functions predict dynamics governed by a 
(3- function pole Xd for rif less than some threshold value (rif < 5 for Z^ 2 ' 1 !; rif < 6 for /^I 1 ! 2 !). 

4.2 Five- Loop Level Results 

It is of interest to examine the stability of the qualitative results described above against five-loop corrections to 
the MS /3-function, corrections which do not enter our determination of ft 2 ^ and /J^ 2 ]. As noted in Section 
2, Pade-coefficients within (3^ 2 \ /3PI 3 ] and /jl 3 ! 1 ! are all seen to be linear in the five-loop /3-function correction f3± 
{(3a/[3o = Ra)- These coefficients, as defined by equations (2.14a), (2.15a), and (2.16a), are respectively tabulated in 
Tables |[ || and ^. In Table |?| we have tabulated the domain of i?4 for which the first positive numerator zero of the 
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n f 


ai 




&2 




Xd 





-5.96723 


-8.28541 


11.0906 


0.16758 


0.15136 


1 


-6.14941 


-8.31070 


10.9765 


0.16262 


0.15007 


2 


-6.68998 


-8.67274 


11.4200 


0.14948 


0.14177 


3 


-8.17337 


-9.95114 


13.2199 


0.12235 


0.11944 


4 


-13.8032 


-15.3432 


20.5809 


0.072447 


0.072160 


5 


70.6188 


69.3580 


-88.9261 




0.79411 


6 


5.95098 


5.02241 


A O ^TO A f\ 

-4.37349 




1.32141 


7 


1.93401 


1.40769 


1.56704 






8 


0.274983 


0.245571 


4.66047 






9 


-0.921639 


-0.321639 


7.31327 


1.0850 




10 


-2.13228 


-0.709208 


10.0353 


0.46898 




11 


-3.60614 


-1.06069 


12.9645 


0.27730 




12 


-5.60030 


-1.43363 


16.1133 


0.17856 




13 


-8.56349 


-1.84921 


19.4406 


0.11677 




14 


-13.6105 


-2.31052 


22.8821 


0.073472 




15 


-24.8114 


-2.81137 


26.3685 


0.040304 




16 


-78.8413 


-3.34126 


29.8352 


0.012684 





Table 6: The coefficients a\, b\ and 62 in ft 1 ^ as given in Eq. (2.13) for n/ = — 16. The column Xd tabulates the 
first positive denominator zero. The column x n tabulates the numerator zero, if positive. For those values of nf 
for which real positive denominator (numerator) zeros do not exist, the corresponding entries for Xd (x n ) are left 
blank. The positive denominator zero Xd occurs before any positive numerator zeros for n/ < 6, suggesting Fig. 2 
type dynamics. When rif > 9, x n corresponds to an infrared-stable fixed point that decreases as nf increases 



Pade-approximant /3-function precedes any positive denominator zero. Such a numerator zero implies a couplant 
with Figure 1 type dynamics, in which the numerator zero is an infrared-stable fixed point. We see from Table |7] 
that such dynamics do not occur at all regardless of R4 unless n/ > 6. An infrared-stable fixed point cannot occur 
for /?[ 2 I 2 1 until n/ = 7 (and then only for R4 < 0), nor can it occur for ft 1 ^ until nj = 9. As rif increases, the 
domain of i?4 for which Figure 1 type dynamics become possible is seen to broaden for all three Pade-approximant 
/3-functions. The overall picture that emerges is quite similar to that anticipated from Tables 5 and 6. In every 
case, dynamics governed by an infrared-stable fixed point (Fig. 1) do not occur below a threshold value of nf, a 
threshold at or above nj = 6. 

Table 8 tabulates the range of R4 for which a positive pole of (3^ 2 \ ft 1 ^ and ft 3 ^ exists and precedes any 
positive numerator zeros, corresponding to the dynamics schematically presented in Figure 2. For n/ < 5, /?[ 2 ' 2 ] 
and /jl 1 ! 3 ! are seen to exhibit such dynamics regardless of the magnitude of R4, the five-loop contribution to the 
/3-function. /3I 3 ! 1 ] exhibits such Figure 2 type dynamics only if i?4 is positive, as the denominator 1 + b\x of eq. 
(2.16a) has a positive zero only if R4 > [b\ = —R4/R3}. Such dynamics, however, become impossible for /3W 3 ] 
and unlikely for f3™ an d (3™ once nf gets sufficiently large, as is apparent in Table 8 from the steadily increasing 
lower bound on i?4 for such dynamics to occur within /3I 2 ' 2 ' and Z^ 3 ' 1 !. 

This large nf behaviour is illustrated for rif — 13 by plots of the independence of the first positive numerator 
and denominator zero 

of £[2|2] ( Fig _ f][i\3] (pig. 7) and/?! 3 ' 1 ! (Fig. 8). The infrared-stable fixed point associated 
with the numerator zero in /?[ 2 I 2 1 and P^ 1 ^ is seen to be less than 1/4 (the assumed threshold for chiral-symmetry 
breakdown pq| ) over the full domain of R4 indicated in Table |^ for infrared-stable fixed point dynamics, a result 
is clearly suggestive of nf — 13 being within the conformal window of QCD. The numerator zero in (3™ is also 
below 1/4, as evident from Fig. 7, until the immediate neighbourhood of its singularity at R4 — 6394. 

The decrease of the infrared fixed point with increasing nf is also common to all Pade-approximant /3-functions. 
This behaviour, as already seen in Tables 5 and 6 for /jl 2 ! 1 ! and P^ 2 \ is illustrated for /jl 2 ! 2 ! in Fig. 9, in which 
the magnitudes of x n are displayed as functions of R4 for nf = {9 — 16}. Such results are consistent with the 
phase structure anticipated in ref. JLlJ , as already noted. It is also worth mentioning that the general picture 
we obtain, particularly the need for a critical number of flavours for an infrared-stable fixed point to occur at all, 
agrees surprisingly well with a lattice study [ p7| . 
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[2|2] 



[1|3] 



[3|1] 




1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



No Ra 
No R 4 
No Ra 
No i? 4 
No i? 4 
No Ra 
No ,R 4 
#4 < 0.21 
Ra < 21.6 
Ra < 57.2 



i? 4 < 62.3 
Ra < 174.9 
#4 < 508.5 
Ra < 1644 
Ra < 6394 
i? 4 < 34187 
Ra < 355521 
i? 4 < 37198800 



Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 



Ra < -105 
Ra < -2.1 



Ra < 1.2 
Ra < 3.6 
i? 4 < 12.9 
Ra < 38.4 
i? 4 < 110 



Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 
Noi? 4 



Ra < 128.4 
Ra < 274.6 
Ra < 594.9 



Ra < 1384 
Ra < 3750 



J? 4 < 340.3 
Ra < 1226 
i? 4 < 6057 
Ra < 92883 



Ra < 14262 
i? 4 < 174589 



Table 7: The domain of i? 4 for which an infrared fixed point occurs for f3^ 2 \ /3W 3 ] and Z^ 3 ! 1 !. The ranges of R 4 
listed (for a given choice of rif) are those for which a positive numerator zero exists and precedes any positive 
denominator zeros of the Pade-approximant /3-function. 



5 QCD's Infrared Boundary 

The case of three flavours is of obvious interest, as Pade-extrapolations to the infrared region can be compared to 
the known empirical dynamics at the onset of the infrared region. We know, for example, that evolution of the 
running coupling constant from its well-determined value at /i = M z pq| leads to a prediction H 



We also know that QCD as a theory of quarks and gluons ceases to exist at momentum scales approaching Aqcd, 
although the interpretation of Aqcd is subject to redefinition for each successive order of perturbation theory. 

Tables 5-7 show quite clearly that an infrared-stable fixed point for rif = 3 QCD is unsupported by all Pade- 
approximant /3-functions considered here, regardless of the magnitude of Ra. This result contradicts the infrared- 
stable fixed point obtained from the analysis of a lower-order expression for the /3-function in refs. |ll| and 
although the absence of such a fixed point at rif — 3 is supported by more recent work jL7) . We also note that all 
Pade- approximant /3-functions considered here exhibit Figure 2 type dynamics, regardless of Ra, except for /jl 3 ! 1 ! 
when i?4 < 0. In such dynamics, the /3-function pole Xd occurs at momentum-scale fi c . As evident from Fig. 2, the 
infrared region fi < \x c is inaccessible to the (real) couplant suggesting that /z c fulfills the infrared cutoff role 

commonly ascribed to Aqcd, the "Landau pole" obtained through use of the truncated /3-function series. 

In Figs. 10-12, we have exhibited the R 4 dependence of the (3- function pole Xd for P^ 2 \ /3W 3 ! and /3' 3 ' 1 ! when 
rif =3. The first positive numerator zero is also displayed in all three figures, and is seen to be larger than Xd 
for all i?4 values considered. We note from Figs. 10 and 11 the apparent stability of Xd in /3I 2 ' 2 ' and fi 1 ^ against 
changes in Ra when Ra is negative. Both figures indicate an infrared-attractor near Xd = 0.4 (a s = 1.3), a value 
well-above the anticipated threshold for chiral-symmetry breaking. 

Given knowledge of an initial value, one can utilize Pade approximant /3-functions to estimate the infrared cutoff 
fi c . To demonstrate this, we assume from the central value of (5.1) that x{\ GeV) — a s (l GeV)/n = 0.153. The 
equation 



a s (n f = 3; \i = 1 GeV) = 0.48 



+0.09 
-0.07 



(5.1) 
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nf [2|2] [1\3\ [3|1] 



All R 4 

1 All R 4 

2 All R 4 

3 All R 4 

4 All R 4 

5 All i? 4 



All R 4 
All i? 4 
All i? 4 
All R 4 
All i? 4 
All R 4 



R 4 > 
i? 4 > 
i? 4 > 
i? 4 > 
i? 4 > 
i? 4 > 
i? 4 > 
R 4 > 
i? 4 > 1.2 
i? 4 > 3.6 
i? 4 > 12.9 
R 4 > 38.4 
R 4 > 110 
i? 4 > 340.3 
R 4 > 1226 
i? 4 > 6057 
R 4 > 92883 



6 R 4 > -105 



i? 4 > -35.4 



7 No i? 4 

8 No R 4 

9 i? 4 > 4848 

10 R 4 > 2964 

11 R 4 > 2990 

12 i? 4 > 3652 

13 R 4 > 5020 

14 i? 4 > 7759 



R 4 > 0.2 
R 4 > 21.5 
No fl 4 
No R 4 
No ,R 4 
No R 4 
No i? 4 
No i? 4 
No R 4 
No fl 4 



15 R 4 > 14491 

16 i? 4 > 174589 



Table 8: The domain of R 4 for which a positive denominator zero exists and precedes any positive numerator zeros 
of /?[ 2 ' 2 1, /jl 1 ! 3 ! and respectively. Such a condition corresponds to the Figure 2 scenario for coupling constant 

evolution, in which the positive denominator zero is an infrared attractor for both a strong and weak phase of the 
couplant. 



attractor Xd = x(p c ): 



We have utilized (5.3) to plot predicted values of \i c for [2|2]; [1 13] ; and [3|l]-approximant /3-functions against the 
unknown 5-loop /3-function coefficient R 4 (= I3 4 /I3q). Fig. 13 utilizes (3^ 2 \ as determined by the n/ = 3 row of 
Table 2, to predict /j, c , given x(l) = 0.153. The curve terminates with fj, c = 1 GeV at R4 — 128, since Xd (the 
infrared attractor) is itself equal to 0.153 at this value of R 4 . What is noteworthy, however, is the stability of fi c 
over the entire range of negative R 4 . Fig. 13 is clearly indicative of QCD's infrared cutoff (mass gap) occurring 
not much below the p mass: fi c — > 660 MeV as R 4 — > —00. 

Thus Fig. 13 is indicative of a lower bound for fi c well-above the phcnomenological value for Aqcd when 
nf = 3. The bound on [i c obtained from fi 2 ^ remains well above Aqcd even if the estimate for a s (l GeV) is 
reduced to the floor of its empirical range. Fig. 14 utilizes /3l 2 l 2 l in conjunction with the lower-bound value of (5.1) 
for a s (l GeV) [x(l) = a s (l GeV)/n = 0.1305]. The figure continues to predict insensitivity to R 4 over the entire 
negative range of R 4l with a somewhat diminished lower bound on /i c : /i c — + 550 MeV from above as R 4 — ► — 00. 

The behaviour described above is virtually identical to that obtained by utilizing /3W 3 ] within (5.3). Figures 13 
and 14 display [i c as a function of R 4 , with x(l) = 0.153 [Fig 15] and x(l) = 0.1305 [Fig 16]. The expression for 
pW3] utilized in the integrand of (5.3) can be extracted from the nj — 3 row of Table 3. Both curves terminate 
at /i c = 1 GeV at R 4 values corresponding to the infrared attractor Xd being equal to x(l GeV) [xd — 0.153 and 
Xd = 0.1305, respectively]. Both curves also demonstrate the same stability of \i c against changes in R 4l as well as 
virtually the same lower bounds for fi c as obtained in Figs. 13 and 14 from ft 2 ^. 

As noted earlier, an infrared attractor associated with an nj — 3 /3-function pole occurs within /jl 3 ! 1 ! only for 
R 4 > 0. In Fig. 17, we plot the /^'^-prediction for the value of fj, c , as obtained from (5.3), against positive values 
of R 4 . Using the central value x(l) = 0.153 from (4.1), we see that fi c > 600MeV over the entire (positive) range 
of Rl. 

For i? 4 < 0, P^^(x) no longer has a pole. This does not mean, however, that a s (/i) has a domain in which /j, 
can get arbitrarily close to zero. Rather, when R 4 < 0, there will exist a Landau pole, i.e. a minimum value of p 
at which a s will diverge: 




(5.3) 




(5.4) 
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Similar dynamics characterize the evolution of a s {n) that follows from the one-loop /3- function \(3{x) = — (3qX 2 ] : 
even though this /3-function itself has neither poles nor non-zero fixed points. We have utilized (5.4) in Fig. 18 to 
find the minimum value of fj, as a function of R4. If x(l GeV) = 0.153, consistent with (5.1), we then find that hl 
approaches 530 MeV from above as R4 — > -co. 

In every case we consider, it is clear that the domain of a s (yu) when rif = 3 is bounded from below by hadronic 
mass scales comparable to or somewhat below the p-mass. Pade-approximant (3- functions appear to decouple the 
infrared region from a s at values of y, substantially larger than Aqc d ■ 

6 Conclusions 

Utilizing Pade-summation QCD /3-functions whose Maclaurin expansions reproduce the known terms of the MS 
/3-function series, we obtain a surprising degree of agreement with infrared properties predicted |fl3| , |l4| , fl7| via the 
't Hooft renormalization scheme |p9f in which the /3-function is truncated subsequent to two-loop order. Within 
the context of N c = 3 QCD, we find clear evidence for a flavour-threshold between rif = 6 and n/ = 9 for any 
possibility at all of infrared dynamics governed by an infrared-stable fixed point. For rif < 6, no approximant-based 
-function (other than the truncated series itself) is able to yield a positive zero that is not preceded by a positive 
pole, regardless of the as-yet-unknown magnitude of the five-loop contribution (R4) to the MS /3-function series 
(2.10). We reiterate that such a zero can be identified as an infrared-stable fixed point only if it is not preceded by 
a positive /3-function pole. 

We also find (Section 4) that when rif exceeds the (approximant-dependent) flavour threshold for possible 
infrared-stable fixed-point dynamics, the magnitude of that fixed point is seen to decrease as n/ increases (Tables 
5 and 6 and Figure 9). The true conformal window of QCD is not expected to begin until the infrared-stable fixed 
point is sufficiently small to preclude chiral-symmetry breaking in the infrared region Jn3ljl4|. In Section 4, we 
corroborate this window's onset at 11-13 flavours, as anticipated from two-loop results JL3|, |14[. 

For values of rif below the (approximant-dependent) threshold for a possible infrared-stable fixed point, Pade- 
summations of the MS /3-function are indicative of infrared dynamics governed by a /3-function pole, as has already 
been observed for N=l SQCD in the absence of fundamental-representation matter fields |l0| . Such a positive pole 
preceding all positive zeros of the /3-function is found to occur for n/ < 5 for all Pade-approximant /3-functions 
constructed from the MS /3-function series, regardless of the unknown five-loop term (R4) in that series, except for 
the /3-function incorporating both a [3|l]-approximant to higher-loop effects and a negative value of R4. Infrared 
dynamics governed by such a pole have a number of phenomenologically interesting properties. Salient among these 
is the occurrence of an infrared cut-off /i c on the domain of the QCD couplant a s (/x)/7r (Figure 2). We find (Figs. 
13-17) the magnitude of \i c to be quite stable against changes in i? 4 , and to be bounded from below by values 
larger than Aqcd and comparable to low- lying meson masses (500-700 MeV). Such an infrared boundary occurs 
even for the [3|l]-approximant case with i? 4 < 0, a case for which a positive /3-function pole does not exist. For 
this case the infrared boundary corresponds to a Landau pole at similar hadronic mass scales (Fig. 18). 

Pole-dominated infrared dynamics have also been argued jl(| to imply the existence of a strong phase that 
devolves to the same infrared attractor as the asymptotically-free phase of the QCD couplant. In such dynamics, 
both phases may share common infrared properties p0| . As we have noted at the end of Section 3, such dynamics 
may provide a clue as to why lattice results for the glueball spectrum appear to be in agreement with results 
obtained in the strong-coupling, large- N c limit of SU (N c ) via supergravity wave equations within a black hole 
geometry. 
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Appendix: Gluodynamics in the 't Hooft Limit 



If the product of N c and a s is finite and nonzero in the N c — > oo limit, the nt — MS SU(N C ) (3- function in this 
same limit is given by |l9| 

^ 2 t4 = -^T A " [1 + + R^ 2 + ^A 3 (+ i? 4 A 4 ...)] (A.la) 

X = N c a s ( f i)/4n, (A.lb) 

Ri = 34/11, (A.lc) 

R 2 = 2857/198, (A.ld) 

R 3 = 86.04326. (A.le) 

The coefficients subsequent to R3 in (A.la) remain unknown at present. However, the known coefficients in (A.la) 
are sufficient to determine Pade-summation functions incorporating [2|1] and [1|2] approximants via (2.12) and 
(2.13): 

[2|1] = _ 11 [1 -2.87219A- 4.00209A1 
^ y > 3 [1 - 5.96310A] V ; 

Bm iX) = 11 A » [1- 5.40935A] 

h K ' 3 [1 - 8.50026A+ 11.8442A 2 ]' 1 ' 

In both of these approximations, the first positive zero of the denominator [0.1677 in (A. 2) and 0.1483 in (A. 3)] 
precedes the first positive zero of the numerator [0.2595 in (A. 2) and 0.1849 in (A. 3)], excluding Figure 1 type 
infrared-stable fixed point dynamics. Rather, both approximations support Figure 2 type dynamics, in which the 
first positive denominator zero (Xd) serves as an infrared attractor for both a weak and a strong ultraviolet phase. 

Such behaviour suggests that both phases may share common infrared properties, as the running couplants in 
both phases evolve towards the same infrared attractor. This behaviour is corroborated by [1|3] and [2|2] Pade- 
summation /3-functions whose Maclaurin expansions reproduce (A.la) inclusive of an arbitrary five- loop contribution 
R4. For example, we find from comparison of (2.15) to (A.la) that 



/3 [1|3l (A) = -yA 2 



1 + aiA 



(A.4a) 



1 + M + M 2 + & 3 A 3 _ 

ai = 15.8424 - 0.0379166i? 4 , (A.4b) 

bi = 12.7515 - 0.0379166i? 4 , (A.4c) 

b 2 = -53.8429 + 0.117197i? 4 , (A.4d) 

63 = -103.614 + 0.184865i? 4 - (A.4e) 

In Figure 19 we have plotted the first positive numerator and denominator zeros of (A. 4a) as a function of R4, the 
unknown five-loop term. A positive denominator zero (Ad) occurs over the entire range of R4 and precedes the 
numerator zero (A n ) when positive, consistent with Figure 2-type dynamics. For the [2|2] case, we find from (2.16) 
and (A.l) that 



/3 [2|2] (A) = _H A 2 

3 



1 + aiA + a 2 A 2 



(A.5a) 



l + M + M 2 . 

01 = 24.5902 - 0.0535237i? 4 (A.5b) 

a 2 = -47.3211 + 0.0844278i? 4 , (A.5c) 

b l = 21.4993 - 0.0535237i? 4 , (A.5d) 

6 2 = -128.203 + 0.249865i? 4 . (A.5e) 

As is evident from Fig. 20, the first positive numerator zero (A n ) is always preceded by a positive denominator 
zero (Xd), which serves as the infrared-attractor for Figure 2-type dynamics, regardless of the unknown five-loop 
term R±. 
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Incorporation of a [3 1 1] approximant within the Pade-summation of (A.l) yields a positive denominator zero 
only if i?4 > 0. This is because = R3/R4, as is evident from (2.16a) and (2.16e). Hence, if i? 4 < 0, infrared 
dynamics along the lines of Fig. 2 are no longer possible. However, we have found that no positive numerator zero 
exists for this regime, also excluding the possibility of Fig. 1 type dynamics governed by an infrared-stable fixed 
point. If i?4 is positive, the denominator zero is once again seen to precede all positive numerator zeros, again 
consistent with Figure 2 type dynamics. All of this behaviour can be extracted from (2.16) and (A.l): 

/3 [3|11 (A) = -Y^ [1 + ai \^fx] +a3X3] (A ' 6a) 

ai = 3.09091 - 0.0116221i? 4 (A.6b) 

a 2 = 14.2929 - 0.0359227i? 4 , (A.6c) 

a 3 = 86.0433 - 0.167698^4, (A.6d) 

61 = -0.0116221.R4- (A.6c) 

The absence of a positive numerator zero when i? 4 is negative necessarily follows from the fact that a\ , a 2 , and a 3 
are all positive when i? 4 < 0. 
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1^ u 

Figure 2: Behaviour of the asymptotically-free running couplant a;(/x) in dynamics governed by a /3-function pole 
at Xd- The point Xd serves as an infrared attractor of both a strong and weak phase of the couplant. An infrared 
cut-off n c necessarily occurs, corresponding to the /3-function pole at x(/i c ) = Xd- 
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Fig. 3 



0.1 




Figure 3: The i?4-dependence of the first positive numerator zero (x n ) and the first positive denominator zero 
(xd) of the [2 1 2] Pade-summation of the QCD /3-function when n/ = 0. The independent variable R4 = (3^/ '/3 is 
proportional to the presently unknown 5-loop contribution (/3 4 ) to the (3- function. The figure shows that Xd and 
x n both exist over the entire range of R4. However, Xd always precedes x n , indicative of the dynamics of Figure 2 
for the evolution of the couplant x(p) from the asymptotically- free ultraviolet region. 
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Fig- 

ii 
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Figure 4: The independence [R4 = /?4//?o] of the first positive numerator zero (x n ) and the first positive 
denominator zero (xd) of the [1|3] Pade-summation of the QCD /3-function when n/ = 0. The figure shows that Xd 
exists over the entire range of i?4- Moreover, Xd precedes x n over the range of R4 for which a positive numerator zero 
exists, leading once again to the dynamics of Figure 2 for the evolution of the couplant from the asymptotically-free 
ultraviolet region. 
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Figure 5: The /^-dependence [i?4 = P4/P0] of the first positive numerator zero (x n ) and the first positive 
denominator zero (xd) of the [3|1] Pade-summation of the QCD /3-function when rif — 0. As in Figures 3 and 4, 
Xd exists over the entire range of R4, and x<i precedes x n where a positive numerator zero exists. 
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Figure 6: The behaviour over a very large region of i?4 of the first positive numerator zero (x n ) and the first 
positive denominator zero (xd) of the [2|2] Pade-summation of the QCD /3-function when n/ = 13. The figure 
shows that x n precedes Xd when i? 4 < 1383, indicative of dynamics governed by an infrared-stable fixed point 
(Figure 1) for this region. The figure also shows that Xd precedes x n when i?4 is very large (i?4 > 5020), indicative 
of dynamics governed by a (5- function pole (Figure 2). The absence of a positive (5- function zero for R4 between 
1383 and 5020 precludes the possibility of dynamics governed by an infrared-stable fixed point for this range of Ra- 
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Figure 7: The i?4-dependence of the first positive numerator zero (x n ) and the first positive denominator zero 
(xd) of the [1|3] Pade-summation of the QCD /3-function when nj — 13. The numerator zero is seen to precede 
the denominator zero over the entire range of i?4 exhibited, consistent with the numerator zero serving as an 
infrared-stable fixed point. This numerator zero ceases to be positive when R4 > 6394. 
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Figure 8: The i?4-dependence of the first positive numerator zero (x n ) and the first positive denominator zero 
(xd) of the [3 1 1] Pade-summation of the QCD /3-function when n/ = 13. The numerator zero is seen to precede the 
denominator zero when R4 < 340, corresponding to the domain of R4 for dynamics governed by the infrared-stable 
fixed point x n . The denominator zero precedes any positive numerator zeros that occur for i? 4 > 340, corresponding 
to Figure 2 dynamics governed by an infrared attractor at Xd- 
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Figure 9: The /independence of the infrared fixed point x n , as indicated by successive plots of x n versus R4 for 
/jl 2 ! 2 ! with values of rif ranging from 9-16. 
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Figure 10: The i?4-dependence of the first positive numerator zero (x n ) and the first positive denominator zero 
(xd) of the [2 1 2] Pade-summation of the QCD /3-function when n/ = 3. The figure shows that Xd always precedes 
x n , indicative of infrared-attractor dynamics (Figure 2). 
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Figure 11: The i?4-dependence of the first positive numerator zero (x n ) and the first positive denominator zero 
(xd) of the [1|3] Pade-summation of the QCD /3-function when n/ = 3. The figure shows that Xd exists over the 
entire range of R4, and precedes x n over the range where such a positive zero exists. Hence Xd serves as an infrared 
attractor. 
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Figure 12: The i?4-dependence of the first positive numerator zero (x n ) and the first positive denominator zero 
(xd) of the [3 1 1] Pade-summation of the QCD /3-function when n/ = 3. The figure shows that Xd exists over the 
entire range of R 4 , and precedes x n over the range where such a positive zero exists, as in Fig. 11. Hence Xd serves 
as an infrared attractor for this case as well. 
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Figure 13: The ^-dependence of /j, c , the minimum value of fi accessible to the couplant x(/j,) derived from /?[ 2 I 2 1 
when rif = 3. The 1 GeV value of the couplant is assumed to be x(l) = 0.153, as discussed in the text 
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Figure 14: The i?4-dependence of /i c , as in Figure 13, except that the 1 GeV value of the couplant is given a 
lower-bound value x(l) = 0.1305. 
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Figure 15: The i?4-dependence of fx c , the minimum value of /U accessible to the couplant x(fj,) derived from /jW 3 ! 
when rif = 3 and x(l) = 0.153. 
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Figure 16: The i?4-dependence of /j, c , as in Figure 15, except that the 1 GeV value of the couplant is given a 
lower-bound value x(l) = 0.1305. 
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Figure 18: The i?4-depcndcnce of the Landau pole /il of the couplant derived from (3^^ when n/ = 3 and x{l) 
= 0.153. This Landau pole (as opposed to an infrared attractor of two ultraviolet phases) occurs only for negative 
values of R4. 
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Figure 19: The dependence on the unknown five-loop term R4 of the first positive numerator zero (A„) and 
the first positive denominator zero (A<j) of the [ 1 1 3] Padc-summation of the SU(N C ) /3-function for the couplant 
A = N c a(n) I '47T in the N c — > 00 limit (n/ = 0). The figure shows that the denominator zero A^ exists over the range 
of i? 4 and precedes the first positive numerator zero A„ when such a zero exists, indicative of Figure 2 dynamics 
with \ d serving as the infrared attractor. 
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Figure 20: The independence of the first positive numerator zero (A„) and the first positive denominator zero 
(Xd) of the [2 1 2] Pade-summation of the SU(N C ) /3-function for the couplant A = N c a(/j,) / Air in the N c — ► oo limit 
(rif — 0). Once again, the denominator zero Ad exists over the entire range of R4 and precedes the first positive 
numerator zero X n . 
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